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Abstract 

Modeling and analysis of cooperative spectrum sensing is an important aspect in cognitive radio systems. In this 
paper, the problem of energy detection (ED) of an unknown signal over Nakagami-m fading is revisited. Specifically, 
an analytical expression for the local probability of detection is derived, while using the approach of ED at the 
individual secondary user (SU), a new fusion rule, based on the likelihood ratio test, is presented. The channels 
between the primary user to SUs and SUs to fusion center are considered to be independent Nakagami-m. The 
proposed fusion rule uses the channel statistics, instead of the instantaneous channel state information, and is based 
on the Neyman-Pearson criteria. Closed-form solutions for the system-level probability of detection and probability 
of false alarm are also derived. Furthermore, a closed-form expression for the optimal number of cooperative SUs, 
needed to minimize the total error rate, is presented. The usefulness of factor graph and sum-product-algorithm 
models for computing likelihoods, is also discussed to highlight its advantage, in terms of computational cost. The 
performance of the proposed schemes have been evaluated both by analysis and simulations. Results show that the 
proposed rules perform well over a wide range of the signal-to-noise ratio. 

Index Terms 

Cognitive radio, spectrum sensing, binary hypothesis testing, likelihood ratio, factor graph, sum-product-algorithm. 

I. Introduction 

It is well-known that most of the licensed spectrum is not fully utilized all the time [1], when fixed spectrum 
allocation is used. Moreover, the rapid deployment of new wireless devices and applications with growing data rates 
creates a spectrum scarcity problem. Cognitive radio networks (CRN) [2] is an emerging solution to the problem 
of inefficient use of allocated licensed spectrum. In this approach, the secondary users (SUs) or cognitive radios 
(CR)s are allowed to sense the spectrum dynamically, identifing the spectrum holes i.e. in the absence of a primary 
user (PU) - in the target spectrum pool and opportunistically utilize it. 
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A. Motivation and Literature 

Spectrum sensing is the first critical step of the CR cycle [2] in order to dynamically utilize the unused spectrum. 
Sensing techniques can be classified as, a) Local Sensing: Each SU individually and/or independently detects 
spectrum holes. Although this kind of sensing is sensitive to fading, shadowing, and model uncertainty, it has a 
simple implementation. A brief survey of different spectrum sensing techniques was presented in [3]. Furthermore, 
it is shown in [4] that the energy detection (ED) is optimal for detecting zero-mean constellation signals, if no 
prior knowledge about PU’s signal is available at the SU, except of the received signal power. Note, that ED is also 
popular due to the simplicity of its implementation [5]. b) Cooperative Sensing: Information from multiple SUs are 
jointly used to detect spectrum holes, and to mitigate multipath, shadowing etc., by exploiting the spatial diversity 
among CRs. It enhances accuracy, reliability, and performance of sensing at the cost of complexity. Moreover, 
cooperative sensing is most effective, when collaborating CRs observe independent fading or shadowing [5]-[ll]. 

Cooperative sensing may be further viewed as distributed detection problem, with the central coordinator to be 
the fusion center (EC). This is also known as centralized cooperative spectrum sensing (CCSS) and is investigated in 
this paper. A detailed survey on the distributed detection was presented in [12]. In this kind of detection, likelihood 
ratio test (LRT) rule is known to be optimal. However, a global optimal solution with coupled local best rules 
is known to be NP-hard and not available in closed-form [13]. Moreover, it was proved that ED is optimal for 
single-sensor detection [14], while identical decision rule is asymptotically optimal for global decision in a large 
network [13]. LRT is implemented using either the Neyman-Pearson (N-P) criterion (maximization of probability 
of detection subject to a constraint on probability of false alarm) or the Bayes criterion (minimization of error) [15]. 

Well-known sub-optimal fusion rules for ideal SU-FC (reporting) channels are: AND, OR, and VOTING [12]. 
Other sub-optimal decision fusion rules over noisy channels are Chair-Varshney fusion for the high signal-to-noise 
ratio (SNR), equal gain combining (EGC) for low SNR, and maximum ratio combining (MRC) for medium SNR 
[16]. Recently, suboptimal fusion rules, known as optimal linear cooperation strategy for additive white Gaussian 
noise (AWGN) and linear-quadratic (LQ) strategy for ideal reporting channels have been discussed in [17] and [18], 
respectively. Furthermore, in [14], [19], the authors used a probabilistic graphical approach to model the optimal 
LRT based fusion for cooperative spectrum sensing. Other sensing algorithms for PU detection, such as optimum 
matched filtering [20], eigenvalue based detection [21], cyclostationary feature detection [22], generalized likelihood 
ratio test (GLRT) based sensing [23], [24], were also reported in literature. 

Good message-passing algorithms, like Pearl’s belief-propagation (BP) algorithm and sum-product algorithm 
(SPA) [25] over suitable graphical models, have been successfully employed for solving inference problems in 
various aplications, e.g. data mining, computational biology, statistical signal processing, and wireless communica¬ 
tions. This approach provides exact solutions for acyclic graphs, while exhibiting a low computational complexity, 
compared to explicit methods [26]. Moreover, BP/SPA is inherently suitable for distributed implementation [25]. 
Therefore, it becomes a practical and powerful tool to solve distributed inference problems, such as cooperative 
spectrum sensing in CRNs [14], [19]. 
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In the spectrum sensing literature, previous studies assume approximate channel statistics [18], [21], [27] or 
known [7], [8], [17] or estimated [23], [24], instantaneous channel state information (CSI). The effects of different 
signal models with known CSI on sensing have also been investigated in [7], [28]. Furthermore, the problem of 
energy detection of an unknown signal over Nakagami-m fading was addressed in a few papers [29], [30], but, 
the results were presented only for high SNRs. Regarding the decision fusion most of the works consider non¬ 
ideal sensing channels with ideal [12], [31] or binary symmetric (BSC)/AWGN reporting channels [17], [19], [32]. 
Optimization of the CCSS scheme over Rayleigh fading and ideal reporting channels, was also addressed in [31]. 
Furthermore, decision fusion over non-ideal reporting channels was introduced by [16], in the context of wireless 
sensor networks (WSN). However, multipath fading on sensing and reporting channels is common in a CRN and 
limits the performance of CCSS. However, none of the above works consider multipath fading on both PU-SUs 
and SUs-FC links, simultaneously. 

Nakagami-m, is a general fading model [33], which often gives the best fit for land and indoor mobile applications 
[34]-[37]. However, cooperative spectrum sensing, in the presence of Nakagami-m fading with channel statistics, is 
relatively less investigated. Moreover, SUs may be mobile in many applications, like object tracking, environment, 
habitat management etc., where the channel estimation is costly. Therefore, spectrum sensing over Nakagami-m 
fading for a wide range of SNR and LRT based decision fusion without knowledge of instantaneous CSI, is useful 
for the system design. Moreover, in a large CRN it also involves conditional and unconditional independence on 
large number of random variables (RVs) and thus it leads to an increase of the overall system complexity. Hence, 
inference over graph with message passing is a good approach for this problem [14], [19]. 

B. Contribution 

In this work, we study the performance of CCSS systems, by assuming that both PU-SU and SU-FC channels 
are Nakagami-m and independent accross the SUs. The LRT statistics is computed through message passing over 
the representative NFG, in order to reduce the computation complexity. Specifically, the main contributions of this 
paper are as follows. 

• Derivation of an LR based fusion rule without knowledge of the instantaneous CSI. Closed-form expressions for 
the system level probabilities of detection, Po, miss. Pm, and false alarm, Pp, are also derived. Furthermore, 
we present an alternate expression for the local probability of energy detection over Nakagami-m fading. 

• Determination of the optimal number of cooperating SUs, needed to minimize the total error rate, Ptot = 
Pm + Pp, as a function of the SNR and the total number of SUs. 

• Modeling of CCSS using NFG and SPA in order to analyse the computation time complexity, compared with 
explicit method. 

C. Structure 

The rest of this paper is organized as follows. Section II refers to NFG and SPA, while Section III represents the 
system model, the assumptions used, and the problem formulation. Expressions for local probalilities of detection and 
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false alarm are presented in Section IV, while the LRT-based fusion rule with NFG-SPA based model, closed-form 
analysis of system-level performance metrics and optimization of the CRN, are presented in Section V. Simulation 
results are reported in Section VI, and the complexity analysis and advantages of NFG-SPA settings are discussed 
in Section VII. Finally, Section VIII concludes the paper and propose some future research directions. 

D. Notations 

Throughout this paper, K denotes the total number of SUs present in the CRN, N denotes the total number of 
(complex) signal samples available for detection, also known as time-bandwidth product. We use to denote the set 
of RVs {oi, ...jO/f}. Here, E[.] denotes statistical expectation, |a| denotes modulus of a, Pa (.) denotes probability 
density function (pdf) and Fa (.) denotes cumulative distribution function (cdf) of a. Nak{m ,.), CAf {.,.), and 
denote Nakagami-m distribution with fading severity parameter to, complex Gaussian, and real Gaussian 
distribution, respectively. 


II. Factor Graph, Sum-Product Algorithm (SPA) 

Probabilistic graphical model (PGM) [38] is an effective way to represent the probabilistic dependencies between 
RVs. Well-known graphical models are Bayesian network (BN), Markov random field (MRF), Tanner graph (TG), 
junction tree (JT), and factor graph (FG) [38]. Among those, FGs are more general, since any BN, MRF or TG can 
be transformed as FG, with no increase in its representation size [25]. Throughout the present work, we consider a 
version of FG, called normal factor graph (NFG) [26], as the PGM. The primary goal of FG-SPA based modeling 
of CCSS is to reduce computational complexity. 


A. Factor Graph 

Factor graph is a standard bipartite graphical representation of a mathematical relation between variables and 
local functions. There are two types of factor graphs [26]: conventional and normal (Fomey-style) factor graph 
(NFG). In an NFG, functions or factors {fj} are represented by nodes and variables {xi} are represented by edges. 
Example: Consider a joint probability mass (density) function /(.) of L variables as f (^xi,X 2 ,X 3 ,X 4 , 
Suppose, the function is factorized as 


f{xi,X2, ■■■,Xl) 


1 


J 

^ {xi,X2,...,Xl} , 


( 1 ) 


J=1 

where Z is a normalization factor. Alternatively, it can be represented through a graph with function nodes and 
variable edges. We consider the factorization with L = 7 and J = 6, where one variable is involved in more than 
two factors. Then, 


/ (a;i, a:2, xs, X4, X5, *6, x?) = ^/a(xi, X2, X3, X4)/b(xi, X5) x 

/c (X2, X7)/o (X4 )/b (xs )/f (xs, X6 ). 


( 2 ) 


Fig.l depicts an example of (2) as a normal factor graph. We summarize the construction of NFG as follows: 
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Fig. 1. Normal factor graph corresponding to Eq. (2). 

• Equality node, 0, indicates variables corresponding to more than two functions i.e. the node with degree, 
D> 2 . 

• Computation of marginal can be performed in an efficient and automated way by using SPA on factor graph. 

• A function can have many factorizations; therefore, it can have many factor graphs. As long as the graphs 
have no cycles, the same marginal will be computed for all. 

B. Sum-Product Algorithm and Message Passing 

Sum-product algorithm (SPA), also known as message-passing or belief-propagation (BP) algorithm, can often 
be applied successfully in situations, where exact solutions to the marginalize product-of-function (MPF) problems 
become computationally intensive [25], [38]. SPA operates over an NFG associated with the global function and 
computes various marginal probabilities by approximating through beliefs. Let us define the message from function 
node fj G Af{xi) to variable edge xi as The message from variable edge xi G ■N'ifj) to function 

node fj is denoted by Mx,^fj{xi), where Af{xi) and Af{fj) are the set of neighboring functions of xi and the set 
of variables involved in function fj, respectively. Message from edge xi to node fj is computed as 

on F[ ( 3 ) 

and message from node fj to edge xi is computed as 

Mf.^,,^{xi)cx ^ fj{Xj) ( 4 ) 

where Af{i)\ {a} denotes all the nodes/edges that are neighbors of edge/node i except for node/edge a. In SPA, 
sum is due to summation and product is due to product operation in (4). In case of continuous variables, summation 
operation is replaced by integration. The proportionality sign in (3) and (4) is used to indicate a normalization 
factor, such that the distribution sums/integrates to one. Final marginal for any variable xi is calculated as belief 
i.e. the product of all incomming mesages as 

b{xi) oc n (xi). 

feMixi) 
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III. System Model for CCSS over Fading Channels 

The block diagram of CCSS system is shown in Fig. 2. It consists of one PU, K number of SUs, and one FC. 
All SUs are independently sensing the PU and then sending their local decisions to FC. Final decision is taken by 
the FC. 



Fig. 2. System block diagram of Centralized Cooperative Spectrum Sensing (CCSS) scheme with K SUs, one PU and one FC. 


A. Assumptions 

Throughout this paper, Nakagami-m fading is used to model rapid fluctuations of the amplitudes of a radio signal, 
and it is assumed that the average sensing duration is much shorter than the average busy-to-idle and idle-to-busy 
state transition periods of PU [31]. Transmit power of PU is assumed to remain constant over a typical sensing 
period and a prior probability of PU’s traffic is unavailable at each SU. 

Furthermore, it is assumed that all SUs stay silent during the sensing interval, such that the spectral power 
remaining in the targeted band is transmitted only by the PU. Next, it is considered that all SUs use same transmit 
power relative to the PU (as in the interweave Dynamic Spectrum Access (DSA) model [39]) and each SU makes a 
binary local decision (hard-sensing) using ED [6]. At FC, the final decision (u) is taken when local decisions from 
all SUs are arrived. We formulate our problem by assuming that all sensing and reporting channels are time-invariant 
(during the sensing process), frequency-flat fading and statistically independent, across different SUs. 

B. Problem Formulation 

Suppose that all SUs monitor the same frequency with the PU (Fig. 2). Spectrum sensing at the fc-th SU can be 
formulated as a binary hypothesis testing problem [15]. The received signal samples at SUk for the two hypotheses 
can be modeled as 


Ho : PU is idle: Xk{n) = Wk{n) 

Hi : PU is busy: Xk{n) = Zk{n) Wk(n), 


(5) 
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where n = k = Wk{n) ~ CAf{0,2a^^) is the sample of AWGN, and Zk{n) represents the 

signal sample, received from the primary user if active. The signal is modeled as RV with average power of 
E[\zk{n)\‘^]=i2z^, which includes the channel gain. In practice, Zk{n) and Wk{n) are independent. 

Classically, the received primary signal samples at SUk are assumed (reasonable approximation for unknown PU 
signals over fading; [14], [23], [24]) to be independent and identically distributed (i.i.d.) complex Gaussian RVs 
with zero mean and variance, This implies that the magnitude of the complex envelope is Rayleigh 

distributed. However, in practice, the received signal at each SUk is composed of a large number of resolved 
multipath components. Therefore, the magnitude of the envelope is the norm of an m-dimensional complex vector, 
where m is the fading parameter of the Nakagami-m distribution [40]. Therefore, \zk{n)\ ~ Nak{m, <Tz^), and the 
average power of the received primary signal can also be written as E[\zk{n)\‘^]=nz^=2ma‘^^ [40]. The average 
received SNR of the PU—SUk link, measured at SUk, is defined as 


(™) _ 


Pp»k 


( 6 ) 


At SUk, the ED computes the energy of the received signal over N samples. The computed energy is compared 
with the threshold r^, which is determined from a given local probability of false alarm, and the binary local 
decision UfeSj—1,1} is generated, where Mfc = —1 and Mfc = l denote absence and presence of the PU, respectively. 
Therefore, the test statistic at SUk becomes 


N 

tk^'^\xk{n)\^ ^ Tk. (7) 

n = l 

Each decision, Uk, is transmitted to the FC over an independent, frequency-flat fading channel. The signal received 
at the FC from SUk is 


Vk = Vk + rik, (8) 

where Uk ^ 2a‘^^) is the observation noise and Vk is the secondary signal, over the SUk—FC link. Similarly, 

we assume \vk\^Nak{m,ay^), where the average power of the signal is =2mcr^^. In practice, Vk 

and Uk are independent and the noise samples at FC and SUs are also independent across different PU—SUk—FC 
links. The average received SNR for SUk — FC link is defined as pU =—j-^. 

Jk 

The vector of received signal at the FC from all SUs is denoted by = {yi,y2, Vk}- For any y^ the binary 
hypothesis problem at FC is 


/o : Primary user is idle: P(j/f'|77o) 

/i : Primary user is busy: P{y^\H\), (9) 

where |/o) and P{yi\Ii) are the distributions of y^ in absence (/q) and presence (/i) of the PU, respectively. 
Final decision (u) is derived at the FC by FRT using these two distributions. 
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We assume that SUs use the spectrum, whenever they detect a spectral hole (white space). The constraint in the 
system model is the probability of erroneous decision about the presence of the PU. Hence, for efficient utilization 
of spectrum, according to Neyman-Pearson criteria [15], the system designer needs to minimize the probability of 
miss or maximize the probability of detection, subject to the constraint that the probability of false alarm satisfies 
a minimum requirement. 


IV. Local Performance Analysis 

This section presents the local performance analysis (probabilities of detection and false alarm at SUk) for the 
system model defined in section-III.B. Throughout the analysis, the variable afc(n) is used as for notational 
simplicity. 

A. Analysis Under the Hypothesis Hq 

Under Hq, the received signal contains only noise and thus, Xk ~ CAf As the number N of complex 
samples is considered as sensing period, the pdf of tk under Hq follows Chi-square distribution with 2N degrees 
of freedom [41]. In this case, the probability of false alarm at SUk is obtained as [29] 

r r(iV, y{N, 

where y{N,x), T{N,x)=T{N) —'j{N,x), and r{N) are the lower incomplete gamma, upper incomplete gamma, 
and complete gamma functions, respectively [41]. Assuming that the noise variance is perfectly known, the local 
threshold Tk is obtained for a given as 




( 11 ) 


B. Analysis Under the Hypothesis Hi 

Under the hypothesis Hi, the received signal contains both PU signal and noise. Therefore, the distribution of 
the test statistic depends on distribution of the envelope of the signal (i.e. |xfe|), which further depends on \zk\- The 
unconditional pdf of \xk\ under Hi is obtained using Bayesian approach, i.e. marginalizing the conditional pdf of 
\xk\ over \zk\- Here, \zk\ is Nakagami-m distributed [33] and is defined as 


T'lzfci(cfc) 


2m—1 / 1 \m <^k 

^Ck / 1 \ 2<t2 

r(m) (,2aiJ ® 


( 12 ) 


The pdf of \xk\ under Hi is obtained as 
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roo 

I (flfc I 1 1 I (fflfc I Cfc ; Hi )P|^^ I (cfc)dcfc 


r(m)iPi ( m, ) + 


r(m)y27ro^ 

aky/2V(m+ |)iPi (m + ^,ak^D^^) 


(13) 


where = 


ZU '-^psu ZU '^psu 

2<t2 “ 25^ ’ ^PSk — 0-2 +^2 — (m) > rj^Sfc 


-('m,') mcr^ 

^ , Ppsl = and iFi(.;.;.) is the 


confluent hyper-geometric function [41, (13.1.2)]. 

As, tk = '^n=i l^kin)]"^ and Xk{n) is complex, tk has 2N degrees of freedom. After a transformation of the 
variables in (13) we get the pdf of tk under Hi as 


PNak{tk\Hl) — 


tk^e 




iPi (m;iV + ftkD^fj. 


(14) 


N\ (2cr2 ' -- -L 

Integrating (14) using Laguerre polynomials [42], [43, (6.9.2.36)], and after some algebra [41, (6.5.12)], the local 
probability of detection can be written as 


(ff) = l_FMak{tk\Hl) 


Pd 


= 1 - 


m y 


2 £7-2 I 


(JV+ 1 ) 


4>2 m, 1; A + 1; 


20-2, ’2cr2„ 


(15) 


A! 

In (15), $2 (■,.;•;•, •) is the hypergeometric function of two variables [41]. The steps for (14)-(15) are given in 
Appendix B. Note that Eq. (15) is more general (as it holds Vm > than [29, (7)], which is restricted to integer m. 
Moreover, Eqs. (15) and [29, (7)] can be implemented via the MATHEMATICA software, which requires truncation 
of infinite series and computation of error bounds [30], [44]. However, a general closed-form expression of is 
intractable, due to the presence of iFi(m;.;.), e“*, and in (14), simultaneously [41]. Therefore, it is interesting 
to find for specific values of m. 


Case-I: m = l 

Eor m = l, (13) is simplified as 


P\x]^\ {klk \ Hl ) ■ 


Aw^. f 27r 


--f 


A,,., 


^k (^( 


(16) 


where = 




and Q{a) is the well-known Gaussian Q-function [41]. Eor the derivation of (16) see 


Appendix C. Setting m = l in (14) and using [41, (6.5.12)] we obtain 
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P{m^l){tk\Hl) = 


(1 ®psfe)< 


(«P 


2<T2,r(Ar) 


-J N. 


ik^p 


20-2^ 


(17) 


where p^s\ = -i^- The cdf is obtained by integrating (17) over tk using [45, (1.2.2.1)]. Hence, the probability of 

'^k 

detection at SUk can be written as 


- / P(m=l){tk\H^) dtu = l- /P(^=l)(4|iyi) 


'ik>Pk 
7 


’ 2ct2 , 


= 1 - 

Case-II: m = 2 

Similarly, for to = 2, (14) can be written as 


r(iv) 


e 2 a^ 
i^PSk) 


7((V,^ 

r(iv) 


(18) 


, N „ 20 

tk e 


P(m = 2) (ffe|77l) = - 


W!(20 


iV+1 


-iPi (2; A1 + 1; . 


(19) 


_(2) 2cr^ 

where pps^=-pr^- With the help of Appendix D, the probability of detection can be expressed as 


Pim=2){tk\Hl) dtu 


= 1 - 






nN-l) r(A-l) 

+ (A—1) (1 — QfpsJ 1 — 


1 + 


Tk 

2A™ 


-2*7 7 ((V, 


Pk^V^h 

2al 

^k 


{apskf m) 


( 20 ) 


Similarly, using other values of to in (14) and integrating over corresponding Pdj,’s can be obtained. For 
example, an expression for pd^., when to= 2, is presented in Appendix F. 

In the next subsection, we present another approach for the determination of pd^. ’s, following the same assumptions 
as in [46], [47]. For the sake of bravity and simplicity, we consider an approximate model, which is valid for moderate 
and high SNRs. 


C. Approximate Complex Representation of the Nakagami-m Envelope 

In [47], it is shown that the exact distributions of real and imaginary parts of the complex signal for Nakagami-m 
envelope are non-Gaussian. As a special case, either the in-phase or quadrature signal may be assumed as zero-mean 
Gaussian, while the other part will be non-Gaussian [47]. In [46], the distributions of in-phase and quadrature parts 
of a signal having Nakagami-m fading envelope, are derived for both to > 1 and 2 < to < 1. However, at high 
SNRs, the proposed model of [46] closely approximates the distributions of real and immaginary parts of signals 
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Fig. 3. ROCs over Nakagami-m fading for different SNRs, with N = 10 and m = l. 


with Rician [46, (59)] for m > 1 and Hoyt [46, (61)], for ^<m<l fading envelopes. These approximations are 
considered here for simplicity. However, for m = 1, Zk{n) can be assumed as CAf (0,2a^^) for all SNRs [14], 
[19], [23], [24]. Therefore, with the help of [46], the complex signals over Nakagami-m fading can be written as 


Zk{n) ~ CAf (0, —^ -f , - < m < 1 

Zk (n) ~ CAf (O, ) , m = 1 

Zfc(n) ~ CAf -f ,m>l, 


( 21 ) 


where = ( "" ^zk cos (ip) and = (yA^) sin (ip), 

6= ‘f defined in [46, (39)]. PU is active under Hi and the distribution of the decision 

statistic is obtained from that of the signal. At high SNR, the distribution of Zfe(n) can be approximated according 
to (21) for different range of m, as in [46]. Thus, the associated probabilities of detection at SUk for different 
values of m are obtained as follows 




N Tfe 

2 ’ 


N Tk 

2 ’ (n.,(l-6)+2a2 J 


Pd 


(m = l) 


j(N, 


\ ^k) 


r(A) 






2Tk 


ilsk + ^cr^k 


( 22 ) 


where Qn {a,b) is the generalized Marcum-Q function [41]. Here, the probability of missed detection is defined 
as pl^k =1 ~Pd^'^ ■ Por '^he derivation of (22) see Appendix G. 
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Fig. 4. ROCs over Nakagami-m fading for different SNRs, with Af = 10 and m = 2. 



Fig. 5. Comparisons of CROC curves over Nakagami-m fading for different SNRs, with = 10 and m = 1, 2. 


D. Comparisons of the Two Models 

Figs. 3 and 4 plot receiver operating characteristic (ROC), i.e. pdk vs pf^. curves for different values of SNR 4 
dB, 10 dB, and 20 dB, over Nakagami-m fading with m = 1, and 2, respectively. 

It is observed from Fig. 3 that both analytical models for m = 1 perfectly match with simulations for all SNRs. 
Further Fig. 4 shows that, for m = 2, at very low SNR none of the models matches with simulations, but, it 
becomes closer to the proposed model in (20). However, at moderate SNRs (4 dB to 10 dB) the simulations match 
with the results of the proposed approaches, while it matches perfectly at high SNR (20 dB). Therefore, we can 
say that for m > 1, the analysis based on the complex signal model [46] is only suitable for high SNRs. However, 
(22) is reasonable and also validates the well-known approximation of Zk{n) for m = 1 [19], [23], [24]. Note that, 
the analytical models proposed in this paper suit well over a wide range of SNR for m > 1. A similar comparative 
study can be presented for | < m < 1. It can also be proved that the proposed model in (61) suits better than the 
complex one in (22) over a wide range of SNR. In this paper, we drop the figure for ^ < m < 1 due to space 


limitations. 
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Fig. 5 plots complementary ROC i.e. vs pf^, curves for different SNRs (0 dB, 10 dB 20 dB) over different 
fading channels (m= 1,2) and compares the proposed model with Digham et. al’s [29]. It is observed that the model 
proposed in this paper is better than [29] over a wide range of SNR. Hence, we always consider the analytical 
models of (18), (20), and (60) for further system-level analysis in the next section. 

V. System-level Pereormance Analysis 

A. Probability Models of the Detection Problem 

The CSS problem, described in Fig. 2, can also be viewed as a system-on-graph (Fig. 6). Thus, we can model 
the problem as inference over the representative NFG and try to solve it by passing messages over the graph using 
SPA. This approach can adopt all the unknown parameters, such as signals, channel effects, noise, and complex 
dependencies in a single framework. It is shown that the FG-SPA approach finds the desired likelihoods in an 
automated way. Accomplishable reduced complexity is shown in Section-VII. 

As explained above, the goal is to find the likelihood functions, P(t/f^|iFo) and P{y^\Hi), in order to compute the 
LRT statistic and thus, to solve the distributed detection problem. These likelihoods can be obtained by marginalizing 
the joint probability distribution of interest over all unknown variables. Thus the detection problem is mapped to a 
Bayesian inference one, in order to find the likelihoods through message passing via SPA over the representative 
NFG. The joint probability distribution P{Hi, uf-, v^, y^) represents the CCSS model of Fig. 2. Likelihood 

functions, represented by P{y^\Hi) for tG{0,1}, can be evaluated as 

P{yi \Hi)=Jp{zi, if, uf, Vi, yi\Hi)dvi duf dt^dzf , (23) 

where P(zf, if, uf, uf, yf |iFi) is the joint distribution of interest. This can be further factorized as 

P(zf,if,uf,uf,j/f|P0 

= P(j/f |uf, uf )P(Mf |if )P(if |Pi, ^f )P(zf )P(uf ) 

K 

= ]^P(j/fc|ufc,Ufc)P(Mfc|ifc)P(ifc|Pi, Zk)P{zk)P{vk)- (24) 

The last line in (24) holds because the channels are independent accross the SUs. Fig. 6 represents the NFG for the 
joint distribution of interest P(2f, if, uf, uf, yf |iTi). Each branch represents the PU-SU-FC path of Fig. 2. The 
equality node, 0, indicates that variable Hi is associated with more than two functions i.e. with all P{tk\Hi, Zk)- 
It computes likelihoods from the joint distribution by employing hard decisions at SUs. The graph has no cycle, 
as all PU-SU-FC channels are statistically independent. Therefore, SPA can compute the exact marginals over the 
graph. 

B. Computation of Messages in NFG-SPA Settings 

In NFG, probability functions are represented by nodes and variables are represented by associated edges. The 
desired likelihoods P(yf |iFo) and P{y^\Hi) have to be computed from the graph. By applying the SPA as message 
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Fig. 6. Forney-style factor graph for joint distribution of interest P(z^, u^, vf^, \Hi)- The graph is shown for two PU-SU-FC channels. 


computation rule [26], intermediate messages are computed and passed between the nodes of the FG. The message 
propagation follows a single step of natural scheduling. As the graph has no cycle, the computation of messages 
starts from the half edge (j/fc) (edge connected to only one node) and leaf node (P{zk) and P{vk)) and proceed from 
node to node. The message from the leaf node (P{zk)) to the connecting edge (Zk) is the marginal 

value of the function (node) with respect to that variable (edge). For half edge, the message (ATj^^_>p(j,^I„j,)) from 
the edge to the node is initialized with the value 1. For intermediate nodes, all incoming messages to that node 
are computed first and every message is computed only once. The messages are indexed with (i),(xi) and are 
shown in Fig. 6 on the corresponding edges of the graph. Dotted arrows show the flow of messages for computing 
the marginals. Step-by-step computations are presented in Appendix A. 

The marginal of Hi on fc—th branch is computed from the final messages of interest similarly as [19, (21)] 

X MHi^p{ti.\Hi,z,,)dzkdtk 

dtkJ Zk 

fTk poo 

= P{yk\uk = -l) P{tk\Hi)dtk+P{yk\uk^l) Pitk\Hi)dtk, (25) 

Jo Jrk 

where 

P{tk\Hi)= f P{tk\Hi,Zk)P{zk)dzk (26) 

J Zk 

is obtained by marginalizing over Zk, and is local threshold at /c-th SU. 

The final message on each branch is computed as the product of messages from node P{tk\Hi, Zk) to edge Hi 
and from edge Hi to node P{tk\Hi, Zk)- The desired likelihood functions for cooperative sensing are obtained from 
the marginal distribution of Hi and computed as the accumulated final message (product over all branches) over 
all edges Hi^^\ Hi^^^ for zG{ 0, 1} from (25) as 

K 

P{yi\Hi) = Y[g{HJ^^). 

fc=i 

Likelihoods are computed in automated way using NFG-SPA. 


( 27 ) 
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C. Analysis of Decision Fusion at FC 

LRT based decision fusion is performed at the FC. It requires likelihoods under FIq and Fli for each received 
signal Uk- Using (25), the likelihood under FIq can be viewed as the message under FIq and can be written as [19, 
(23)] 


P(yk\HQ)^Mp(^tk\zk,Ho)^Ho 

= P{Vk\uk = -l){l - Pff) + P{yk\uk = l)P/i,. ( 28 ) 

Similarly, the likelihood under FIi can be written as 


P{yk\Hi) - 

= P{yk\uk = -1)(1 -Pdfc) + P{yk\uk = l)pdj,. (29) 

As the observations are independent, the LRT for choosing FIi in N-P method can be written with help of (27) as 


r K. Pjyil Hi) _J^ P(yk I -1) (1 -Pd^)+P(yk I l)pdk 

p{y^\Ho) l}p(yk\-m-pfk)+p<ykmpfk ’ 

where A is the threshold at FC. In N-P settings, this is obtained by solving 


(30) 


/ P{yi\Ho) dyf = rP{L\Ho) dL = PF, 

JL{yf)>X 

for a constraint on probability of false alarm at FC. 

According to (8), Uk ~ CAf{0,2a'^^) and the envelope of the received signal at FC over SUk — FC link is 
Nakagami-m distributed i.e. |vfc| ^Nak{m,av^). In general, P{yk\uk) follows (13) by replacing Xk, and crj,^ 
with yk, and a'^^, respectively. For to= 1, as P{yk\uk) follows (16), the LRT statistic (using Appendix E) can 
be written as 


Vk 

9 — 2 

K e + 


{pd^ (SnJ/fc)} y^Bn^ykC 


Vk 

■ 2An 




L{yi) = ll 

'“-1 e ^<k +{pf^-Q{B„^yk)}s/^Bn^yke 

Similarly, for m = 2, the LRT statistic (using Appendix E) is given by 


(31) 


Vk 

2An 


Vk 

2An 


JFB^^yke ^ 

^ ^ Bn^VkC +{pf^-Q(B„yyk)}V^Rnke 


(32) 


_ 




where A - -C rr^ ] - f 1 -C R - ^ 

Where - (cr^^ + a„J - ^ Bnk - , P,/, - ^2^ , a./, - 

i?„j,=(2 + B'^^yk^), and pf^ is obtained from (10). The Pdf& are obtained from (18) and (20) for m = 1 and 2, 

respectively. Similarly, other likelihood ratios may be obtained by substituting corresponding values of m in (30) 

with associated The LRT statistic for m = | is derived in Appendix E. 
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According to this approach, local thresholds are derived, based on the probability of false alarm, and therefore, 
CSI is not needed at each SU. However, local probabilities of detection and final LRT statistics, i.e. L{.) values, 
depend on channel statistics instead of instantaneous CSI. Moreover, system-level thresholds for LRT at FC are 
selected based on the L{.) values. Therefore, we can state that, for the overall system design, the proposed fusion 
rule requires the knowledge of channel statistics, i.e. average SNRs and instead of the instantaneous CSI. 


D. Closed-form Analysis 

Eqs. (31) and (32) are LR based fusion rules, but, the derivation of closed-form expressions for Pd and Pp from 
them, seems to be a hard task. However, LR based optimum fusion rule can be approximated as Chair-Varshney 
one [13], under the assumption of high SNR and identical detectors (i.e. pd^ = pd, Pf^ = Pf V/c) [16]. It is already 
assumed that yfs are statistically i.i.d. for large number of SUs. To derive the closed-form expressions for Pp and 
Pp, we also define as Ki the number of SUs for which j/fe > 0 and K — Ki the number of SUs for which < 0. 
Hence, the log-likelihood-ratio (LLR) from (30) can be written in terms of Ki and K 


log{L{y^)) = {K-Ki)Log 


P{yk<0\ - l){l-pd)-\-P{yk<0\l)pd 


KiLog 


P{yk<Q\ - l){i-pj)+P{yk<Q\l)pj 
Pjyk > 0| - l)(l-pd) -f Pjyk > 0\l)pd 
P{yk > 0| - I){l-Pf ) -f P{yk > 0|l)p/ ■ 


(33) 


Therefore, Ki is binomial {K,p) distributed, where the probability of success p is defined as p = P{yk > 0). 
Let us denote, pi and pq are the probabilities of success under Hi and Hq, respectively. Then, the system-level 
detection performance, i.e. Pp, Pp and Pm{= ^ — Pd) can be computed using the binomial distribution. Associated 
closed-form solutions are 


Pd^J2 ( 34 ) 

Pp = ^ po-’(l(35) 

where Ki takes values from 0 to K. For m = 1, pi and po can be evaluated as (see Appendix E for the derivation) 




{Pd - 

{pf - |)cr„fc 
V^Vk+^k 



(m=l) 

Pd 







(36) 


Similarly, pi and po are obtained for other values of m with associated pd (see Appendix E and E). Hence, we 
state that N-P test is equivalent to J-out-of-K rule when 


• All PU-SU-EC channels are statistically independent. 


• All SU’s are employing identical decision rules and transmitting hard local decisions to EC. 


The theoretical ROC curve at EC is obtained from (34)-(35). 
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Fig. 7. Number of cooperating SUs (/) vs Ptot curve when all channels as Nakagami-1 faded. pj^=0.03, A^ = 20, K= 20 and 15. 

E. Optimization of Cooperative Spectrum Sensing 

In general, when Pp increases, Pp also increases and as a result Pm decreases. However, an effective system 
design should always minimize the total error rate, Ptot = Pm + Pp- An optimal voting rule that minimizes Bayes 
risk function has been addressed in [13, (pp. 94)]. Minimization of Ptot over ideal SU-FC channels using some 
linear approximation has also been reported in [31]. Next, we investigate the optimal number, lopt, of SUs required 
to minimize the total error rate over non-ideal SU-FC channels. 

Let, I be the number of SUs required to achieve certain values of Pm and Pp. From (34) and (35), let us define 
a function D{1) = Pp — Pp,. Therefore, 


D{1) = (Pm + Pf-1) 



(37) 


From the properties of ROC curve [15], it is known that > pf. Then, from (36) and (58) for m = 1 and 
2, respectively, holds that pi > pq. Under this condition, it is evident from (37) that minimization of Ptot is 
equivalent to the minimization of D(l). Now, when I is increased by 1, it holds that 



(38) 


From (38), D(l -f 1) > D(l) if. 


i/1 \K-l ^ i/i \K-l 

Pi (1-Pi) >Po(l-Po) 


or, I X In— >iK-l)x In-—— 
Po 1 - Pi 


K in— 

or, I > -- , where /3= /° . 

- 1 + S ’ ^ 

i-Pi 


(39) 


Similarly, when I is decreased by 1 
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D{l)-D{l-1) 

From (40), D{1) < D{1 - 1) if 


Z-1/1 \K-l + l ^ \K-l+l 

Pi (l-pi) <Po (1-Po) 

or, (i - 1) X In— <(K-l + l)x In ^ ~ 

Po 1 - Pi 

Hi 

or, / < 1 H-- , where B = — , . 

1 + /3 ln\^ 

i-pi 

Thus, from (39) and (41) the optimal number of SUs, which minimizes Ptot can be written as 


(41) 


lopt - 


K 

1 + / 


(42) 


where 


InJ^ 

B= _^ 

i-pi 

and pi, Po are given by (36) or (58). Therefore, lopt depends on K, po, and pi, subject to the constraint pi > pg. 
These probabilities of success under Hi and Ho further depend on the average SNRs of SU-FC channels, local 
thresholds (r^), noise variance (ct^j.) of individual SUk’s., and average SNRs of PU-SU channels. Hence, it is 
observed that the optimal number of SU (lopt) depends on K, r, and average SNRs jjps and p^/. Moreover, it is 
concluded that a fast spectrum sensing algorithm can be executed by considering only lopt SUs instead of all K in 
a CRN, with statistically independent channels. 

Fig. 7 plots the number of cooperating SUs vs Ptot for different values of K and different channel conditions. 
The figure shows that, for a fixed average SNR of 10 dB over SU-FC channels, the optimal number of cooperating 
SUs increases from 11 to 14, while the total number of SUs increases from 15 to 20, respectively. It can also 
be observed that as K increases, Ptot of corresponding lopt decreases. However, to achieve a given Ptot (say, 
0.63) with fixed K = 20, the required number of cooperating SUs (1) increases from 10 to 12, as the average SNR 
decreases from 10 dB to 7 dB. It infers that, for a fixed values of K, the value of I increases to achieve a given 
(bounded) Ptot, as SNR decreases. Hence, we can say that there exist an optimal number of cooperating SUs for 
different K and different SNRs, subject to minimization of Ptot- 


VI. Simulations and Discussions 

In this section, simulation results are presented to evaluate the proposed CCSS scheme. The parameters used 
in the simulations are: K = 5, and 10, TV = 20 and number of Monte Carlo iterations = 10, 000. Independent, 
frequency-flat, Nakagami-m fading with m = 1 and 2, and cr2,^=(T^j,= 1, are considered. In the following, ED 
threshold (r^) at SUk, V/c, are obtained from (11), to maintain a target local probability of false alarm p/^=0.03. 
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Fig. 8. Comparison of ROC curves for various fusion rules over Nakagami-1, 2 sensing and reporting channels where mean of {pps} ~4 dB 
and {psf} ^4 dB with K = W, N = 20, pf = 0.03. 



Fig. 9. ROC curves over Nakagami-1 fading for different means of {psf } ~ 2 dB, 7 dB with fixed mean of {pps} ~ 4 dB, K = 10, and 
Pf = 0.03. 


However, each SUp has different probabilities of detection based on the different valuse of m (e.g. (18), (20) 
for TO = 1,2, respectively). The system-level performance is quantibed by the ROC. The average SNRs of sensing 
and reporting channels are dehned as (dB)=101og]^Q and (dB)=101og]^Q ^. 2 "'° ^ Vfc, respectively. 

Fig. 8 describes the ROC curves at different LRT threshold levels, for different fusion statistics of the CCSS. It 
is assumed that, {-4, -2,0, 2,3, 5,10, 8, 7,11} dB and {ps/{-5, -3,-1,0, 2,4, 7,12,10,14} 

dB, m= 1,2, and K — 10. It can be seen from Fig. 8 that the performance of the LRT based proposed fusion rule 
is better than the two suboptimal schemes, e.g. EGC and MRC [16], when to = 1. Fig. 8 also depicts that the ROC 
curve improves as to increases for all channels, since the effect of fading decreases. 

A more practical heterogeneous scenario with 10 SUs is also considered in Fig. 8. It is assumed that in the first 
six PU-SU-FC channels, to = 1 and for the rest four channels, to = 2. It is shown that ROC performance improves 
for the heterogeneous case, when compared with that of to = 1, because of the SU’s spatial diversity (i.e. different 
fading severity parameters). 


Fig. 9 represents the ROC curves at FC of the LRT based proposed fusion rule, with to = 1, AT = 10, and 
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Fig. 10 . Probability of detection as a function of mean of avg. SNRs. Mean of {psy}?iil0 dB and that of {pps} varies from —20 dB to 30 
dB. 7 ^ = 10 , A^ = 20 , py = 0 . 03 , and Pf = 0.02. 


Pf = 0.03. We consider single set of average SNRs for PU-SU links as {Pps^.} {~4, —2,0, 2,3, 5,10,8, 7,11} 

dB, where the mean of all the average channel SNRs for the PU-SU links is si 4 dB. However, the two sets (low and 
high) of average SNRs for SU-FC links have been considered as ~ {—10, —5, —2, —1,0,3, 5, 7,10,12} 

K 1 U k=l 

dB and Ipsf^. | ~ (—3,—1, 0,4,6,8,10,12,14,20} dB. Here, the means of all the average channel SNRs for 

SU-FC links are Ri 2 dB and 7 dB, respectively. The system-level analytical results and Monte Carlo simulations 
are obtained from (34)-(35) and (31), respectively. As expected, ROC performance increases as average SNRs of 
SU-FC channels increase. 

For better clarity, in Fig. 10, the local and system-level probabilities of detection are plotted as a function of the 
mean of the average SNRs of all the PU-SU channels. The cases of m=l,2 with K=1Q, p/ = 0.03, and Pf = 0.02, 
are considered. Here, SU-FC channels have mean of all the average SNRs r: 10 dB and that for PU-SU links are 
varying from —20 dB to 30 dB. Note that there are significant variations for average PU-SU channel SNRs. In 
Fig. 10, the performances of both cooperative and non-cooperative sensing schemes are presented for different m 
values. It shows that cooperation among SUs significantly improves the probability of detection compared to the 
non-cooperative case, over a wide range of SNR. It is also observed that for m = l, the coopeative scheme achieves 
0.95 probability of detection at 5 dB of SNR, whereas non-cooperative sensing reaches to the same at higher SNR 
(17 dB). Fig. 10 also shows that the probability of detection increases significantly for all SNRs, as m increases 
from 1 to 2, in both the cases. 

Fig. 11 plots ROC curves obtained from (34)-(35) for different values of AT = 5,10 and m=l,2 with pf =0.03. 
The case of {pps}R{—4, —2,0, 2, 3, 5,10, 8, 7,11} dB and {/Os/}«{l, 3, 5, 8, 9,10,13,15,17,9, } dB are considered. 
Here, the mean of the SNRs for PU-SU and SU-FC channels are 4 dB and 9 dB, respectively. For K = 5, SNRs 
are drawn from index 3 to 7 of the above sets (to maintain the same mean for both K). It shows that detection 
performance increases significantly as the number of cooperating users increases from 5 to 10 for m = 1 and 2 
both. The reason is the accumulation of information from more number of SUs. 
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Fig. 11. ROC curves at FC for different K{5, 10) over Nakagami-1, 2 channels with mean of {pps}~4 dB, {psf }~9 dB, N =20, Pf= 0.03. 


VII. Complexity Analysis and Advantage oe NFG-SPA Based Approach for the Sensing Problem 

The usefulness of the NFG-SPA based approach to cooperative spectrum sensing problem can be analyzed via 
the computational cost to find the marginal likelihoods. A function may be factorized in several ways, resulting in 
respective FGs. As long as the graphs are acyclic, the same marginals will be computed. The physical meaning 
of acyclic is that all PU-SU-FC channels are statistically independent. However, in correlated cases it may lead to 
more complex SPA, due to the presence of cycles in the graph [25]. 

The complexity of an algorithm is usually measured in terms of its performance as a function of the size of the 
input. Here, it depends on the number of functions (nodes) and variables (edges) [26]. Let us consider a function with 
M variables that may be factorized (acyclic) in F factors. We assume that each variable is defined over the same 
X discrete-sampled values from a continuous domain or domain X for discrete RVs. The function is represented 
by an NFG of F nodes. Suppose, among F nodes, di have degree 1, ^2 have degree 2,..., du have degree D, 
where D is the maximum degree of a node in the graph. For simplicity, we assume that 1 CPU cycle is required for 
computing message at any node and that time is negligible at any edge. Then the total complexity (in CPU cycles) 
in the graphical method can be written as [26], CpG='^f^iidi\X\^, where \X\ is the cardinality of the domain 
of RVs. Following the conventional (explicit) method, it is required to integrate out all other variables which has 
a complexity equal to Ccn = M\X\^ . Thus, the NFG-SPA based approach (0(|X|^)) is computationally more 
efficient than the explicit method {0{\X\^). In general D « M, hence Cfg«Ccn- 

Example: Consider the graph of Fig. 6 with K=l, i.e. the case of M = 6 and assume |2f| = 2. Therefore, the total 
complexity to find the marginals using the graphical method is 

Cfg= [Ix2x2^+2xlx2^+3x2x2^] = 60 cycles. (43) 

In comparison, the complexity in the numerical method is 

CcN = 6 X 2® = 384 cycles. (44) 

It is interesting to note that NFG is computationally more efficient by ^ Ri6.4 times than the conventional method. 
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Fig. 12. Computation complexity vs. number of RVs involved in computations over Nakagami-1 faded PU-SU and SU-FC channels. 


To understand the advantage of the factor graph for the distributed detection problem addressed in this paper. Fig. 
12 plots the computation complexity as a function of the number of RVs involved. The desired likelihood P{yk\Hi) 
is obtained by marginalizing P{zk,tk,Uk,Vk,yk\Hi) in (23). Each variable of interest of the function is defined 
over the same X discrete-sampled values from a continuous domain (i.e. uniformly quantized with Q= |2f| =4 
levels). We consider two cases as |V| =2, and |2f| =4 for K = 10. From the previous discussions and Fig. 6, in 
this case M = A, D = 2. The total complexities to find likelihoods in both conventional ((23)) and graphical ((24)) 
methods are computed as 

CFG = ^idi\X\'=[1x2x2^+2 x3x2^] x 10=280 cycles ] 

i=i ? , 

Ccn = M\X\^ =[4 X 2*] X 10 = 640 cycles J 

for |X| = 2, and 

D \ 

CtG = '^idi\X\^ = [lx2xA^+2x3x4‘^]xlO=10A0 cycles 

^ , 

Cj^=M|X|*^=[4 X 4^] X 10 = 10240 cycles 

for |X| = 4, respectively. Fig. 12 shows that the complexity monotonically increases with the number of RVs 
involved (i.e. as the system becomes more complex) in both methods. It also shows that complexity increases as the 
quantization level Q increases. However, it is interesting to note that the rate of increment in NFG-SPA setting is 
considerably less than that of the explicit method. It implies that NFG-SPA setting is computationally more efficient 
(e.g. ||^Ri2.4 times for |Ar|=2) than the conventional method (for acyclic FG). Therefore, we can state that for 
a large network and higher quantization level, graphical method is more effective from an implementation point 
of view, even for statistically independent channels. Moreover, collaborative spectrum sensing is very effective, 
when collaborating CRs observe independent fading [5], which results in acyclic FG. Hence, NFG-SPA settings are 
suitable from implementation perspective. 
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VIII. Conclusions 

In this paper the problem of centralized cooperative spectrum sensing over Nakagami-m fading in a CRN, was 
addressed. We have presented a new fusion rule based on LRT, which requires only the statistical characteristics of 
the wireless channels between PU-SU-FC. It was compared with other suboptimal rules e.g. EGC, MRC and better 
performance has been observed. Furthermore, we derived a novel expression for the local probability of detection 
for ED over Nakagami-m fading. It is shown that the proposed models perform better, over a wide range of SNR 
for different values of m, compared to the approximate complex signal representations. Closed-form solutions for 
the system-level probabilities of detection (Pd) and false alarm (Pp) were derived. Eurthermore, an expression 
for the optimal number of cooperative SUs, needed to minimize the total error rate, is also obtained. It leads to a 
fast spectrum sensing algorithm, by considering only the optimal number (lopt) of cooperative SUs instead of the 
total one. Eurthermore, decision fusion based on SPA over representative factor graphs was used. It was shown that 
a considerable amount of gain in computational complexity can be achieved through NEG-SPA settings, even if 
PU-SU-EC channels are assumed to be statistically independent. More complex channel conditions, e.g. correlated 
PU-SU-EC channels, and design of fusion rules for soft decisions can be investigated as future works. 


Appendix A 

Message Passing on Eactor Graph of Eig. 6 


— 1 ; — P{vk) 

= P{Vk)', 

(irj)/ P(2/fc|rife, u/j;)P(pfc) dvk P(yfcl'Ufc) 

^P{yk\ - l)I{tk <Tk)+ P{yk\l)I{tk> Tk) 
{vii)Mti^^P^ti^lHi,zi,)=P{yk\ - l)I{tk<Tk)+P{yk\l)I{tk>Tk) 

= P(zk) 


{ix)M^^^Pi^tk\Hi,zk) = P{zk) and {x)MHi^p{t^\Hi,zi,) = 1 
{xi)Mp(ti,\Hi,zk)^Hi = / Pitk\Hi,Zk)dzkdtk x (ix) x (t 

= 11 Pizk)P{tk\Hi,Zk) [P{yk\uk = -l)/(4 < Tk) 

J ti.J Zl. 


P P{yk\uk = l)I{tk > Tk)] dzkdtk 


Appendix B 

Steps from (14) to (15) 

Integrating (14) using Laguerre polynomials [43, (6.9.2.36)], [42], the cdf of tk is obtained as 


FjVafc {tk\Hl)= M 

1=0 


(QpsJ'(m). 

i\ iN + i)\ ^ 


N+l+i, 



(45) 


(46) 
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In (46), {G)i = for i = 0, l,2,..oo is the Pochhammer symbol [41]. Using [41, (6.5.12)] in (46) and then 

[41, (13.1.2)], we get 


where M = 




^Nak (ik |ffl) — 


(AT + i)! 

i=0 


(JV+l-H) 

iFi (l;At+l + i + l;;^ 


(At + 1 + i) 


M 


2 <t 2 „ 


(iV+1) 


At! e 


<1>2 m, 1; At + 1; 


Oipsf^'^k Tk 


2al. ’2al 


and (a)! (a+1)^- = At! (At + l)j_|_j, when a = N + 1 + i. 


(47) 


Appendix C 
Derivation of Eq. (16) 

For m=l, using the values of iFi |,a), iFi (l, |,a), E (|) [41] , and after some algebra, (13) may be 
simplified as (without loss of generality, index k is dropped from variables A,B,D) 


Similarly, for m 
becomes 


d(1) 




AniV^ 

2 


ae“ ( — 


i, V Au,, 


+ 1 


- 


e Q {—B^a). 


(48) 


Ajuy/^ A, 

2, using the values of E (|), iFi (|, |,a), and iFi (2, f,a) [41] , and after some algebra, (13) 


d(2) 

\^k\ 


(a|J7i) 


^in Bin a 


aBinC 


2 A ,„2 I ^ 

where and + tr^,J . 


- + (2+77,„a^e — BmOi) 


(49) 


Appendix D 
Derivation of Eq. (20) 

As 'y {N,a)=N~^a^e~°'iFi (1; At + l;a) [41, (6.5.12)], and using [41, (13.4.3), (6.5.12)], we can write 
lEi (2; At + 1; a) = (1 - A)iEi (1; At + 1; a) + NiFi (1; A; a) 

= ^^ e°{a 7 (A-l,a)- 7 (A,a)} . (50) 

Therefore, (19) can be simplified as 

P(m=2) {tk\Hl) = 

_ {^~°‘P‘k)*k 

{l-a,nfe~ 

2ai^r{N-l){a,s,f 

Itegrating (51) over tk by parts, then, using J 7 {N, a) da = a'j {N, a)+E (N + 1, a) and [45, (1.2.2.1)] we get the 
local probability of detection for m = 2 in ( 20 ). 


'vsk'^k^ nv—1 


2^ik 


‘2<ylk 


-7 A, 




‘2<^lk 


(51) 
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Appendix E 
Derivation of Eq. (35) 

Consider, for all SUs and are same i.e. pd*, =Pd and pf^, =pf Mk. As G {—1,1}, P{yk\uk) follows 
(13). Thus, for m = 1 P{yk\uk) follows (16), by replacing XkT(T'l^ and cr^,^ with yk,cry^ and cr^^, respectively. 
Then, 


Uk^nf^PriVk 


E(pfc|Mfc)= ~'y + 

AW2tv 


An 


e (^( UkPnyk^ • 

V 27r 

Therefore, likelihood under Hi is written as 

P(m^i){yk\Hi)= P{yk\uk = -1)(1 -Pd) + P{yk\uk = l)Pd 

e -\-{p^ — Q{Bnyt^}y/^Bnyke~^^ 


A„v^| 

where A„=(cr^^+ ct^^) and Bn= ^ further compute 


poo 

Pi^P{yk > 0|ffi)= / P{yk\Hi) dyk 
Jo 


An^PpK LV 2 

(m=l) 


^cr„j, +PdV^f3nA„ — 


V^Bn 


(An (Tyj^'yAn. 


(ri_ 

2 


2/^'^fc 


Eor m = 1, by replacing pd with p/ in (53) and (54) we get P(m=i){Vk\HQ) and po, respectively as 


Po = ^(yfc > 0|i^o) = 7 ,+ 


1 , (P/ - i) 




2 


Similarly, for m = 2, P(pfe|Mfc) follows (49). Thus, we get 


(52) 


(53) 


(54) 


(55) 


P{yk\uk) = 


BnPk 

2Al 


^ +ufc (2 + £-2^^ Q i-UkBnyk) 

V ZTT 


and the likelihood under Hi as 


P(m=2){yk\Hl) — 


^nPnyk 

2V^Al 


JLk 


BnPkC ^'^’'k+{pdi^-Q{Bnyk)}V^ {2+Bnyk) e ^ 


yk 


Then, pi is computed from (57), using [45, (1.5.3.1,8)], as 


Pi- 




+ (^ni, + 


b! 


(m=2) 1 5 




(56) 


(57) 


( 58 ) 


Similarly, for m = 2, by replacing pd with py in (57) and (58), we get P(jn= 2 ){yk\HQ) and po, respectively. 
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L{yi)^Y[ 


(m=l) 


1-1 


-^erf 






1 + V -^erf 


y2<(i+2p./,)yj 


fc=i (1-p/J 


1“ V -^er/ 




+ P/)c 


1 + V -^er/ 


\/2-».(l+2+/.) 


(62) 


Appendix F 

Derivation and the LR eor m= \. 

As discussed earlier, E\\zk^]=2m(j1^ and i<^[|wfc|^]=2mcr^^. Setting m = | in (14), we get 


+ 

tk e 




2+A 


(4) 


(A)! (2o 


>iV+l 


X iFi [ A + 1, tkDn,^ 


(59) 


-(-) 

where ppi' = is the average SNR of PU — SUk link. Now, integrating (59) over tk and using [41, (6.5.12)], 
the probability of detection for m = 4 may be computed as 




r(A) 


(60) 


As, UfcG{—1,1}, P{yk\uk) is obtained by replacing Xk, crj, , and cr^ with yk, , and , respectively in (13) 


P{yk\uk) = 




+ <Xn 




1 + A/^er/ 


ykUk 


\/2cr4(o-2^ +cr4) 


(61) 


where A„ = (cr^^ + (1 + ) and average SNR of SUk — FC link is defined as ^ Now, 

LRT statistic may be written as (62). For both m = ^ and 2, LRT thresholds (A) are obtained by (numerically) 
solving P(Hffo) dL = PF, for given probabilities of false alarm. 


Appendix G 

Derivation oe Eq. (21) and Eq. (22) 

Case-I: ^ < m < 1: In this case, the complex Nakagami-m envelope are expressed as Hoyt approximation 
[46, (61)]. The real and imaginary parts are Af (^0, Af ^0, ^, respectively, where 

Therefore, Zk{n) ~ CM ^0, , Now, as Wk{n) ~ CM (0,we may write 

Xk{n) ~ CM ^0, F j. The test statistic tk is the sum of N squares of i.i.d. M ^0, 

and another N squares of i.i.d. M ^0, ^+c2,fc j j real Gaussian RVs with each having zero mean. Thus, the 

pdf of tk under PIi is the sum of two central Chi-square distributions with each having N degrees of freedom 
(DOE) and written as P(^i<:m <i)(ffc|T^i) — 


r (f) (Hz, (! + &) + 20 ^ r (f) O (!-?>) + 20 ^ ' 
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Hence, for ^ < m < 1, the corresponding probability of detection at each SUk is obtained by integrating (63) as 


Pd. 


7 if 


’ (n.,(i+6)+2a2 


2 ’ (n,^(l_6)+2a2 J 


r(f) r(f. 


(64) 


Case-II: m = 1 : This is the case of Rayleigh distribution where Zk{n) ~ CAf (0,20-^^,) [28]. Then, the signal is 
written as Xk(n) ~ CAf (0,2 + 0-2,^.)). Therefore, the pdf of test statistic under Hi follows central Chi-square 

distribution with 2N DOF. The probability of detection at SUk is computed by integrating P(rn=i){tk\Hi) as [19] 


(m=l) H 

Pdk =!-■ 


7 AT, 




r(iV) 


(65) 


Case-Ill: m> 1 : In this case, the complex Nakagami-m envelope are expressed as Rician approximation [46, 


(59)]. Here, the real and imaginary parts are Af j and Af distributed, respectively. Therefore, 

Zk{n) ~ CM j/iQ,^,H.j,^.Now, as Wk{n) ~ CAf (O, 2cr^J, we may write Xkin) ~ CAf + 2cr^^ 

Then, the test statistic (tk) is a sum of 2N squares of independent and non-identically ditributed Gaussian RVs with 
each having non-zero mean. Therefore, the pdf of tk under Hi follows non-central Chi-square distribution with 


2N degrees of freedom and non-centrality parameter = l]f=i 
With help of [41], it may be written as 






^Qz 




+ 1 


P{m>l)itk\Hl) = - 


(S;) 


N-1 i __^+U2 

(n,J 2 


-fiV-l 


Hzi,2tk j 

I 


( 66 ) 


where 0<ffe<oo, =^sk + and /^(.) is the modified Bessel function of the first kind of order iz. The 

probability of detection at SUk is given as [28, (2.13)] 


pliT"^ = QN\^k, 


2Tk 


Qs. -I- 20-2 


(67) 
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